MODELING USUAL AND UNUSUAL 
ANISOTROPIC SPHERES 



Stefano Viaggiu 
Dipartimento di Matematica, Universita "Tor Vergata" 
Via della Ricerca Scientifica, 1 

Rome, Italy 00133, 
viaggiu® axp . mat . uniroma2 . it 

April 6, 2009 



Abstract 

In this paper, we study anisotropic spheres built from known static 
spherical solutions. In particular, wc arc interested in the physical 
consequences of a " small" departure from a physically sensible config- 
uration. The obtained solutions smoothly depend on free parameters. 
By setting these parameters to zero, the starting seed solution is re- 
gained. We apply our procedure in detail by taking as seed solutions 
the Florides metrics, and the Tolman IV solution. We show that the 
chosen Tolman IV, and also Heint Ha Durg IV, V perfect fluid solu- 
tions, can be used to generate a class of parametric solutions where 
the anisotropic factor has features recalling boson stars. This is an in- 
dication that boson stars could emerge by "perturbing" appropriately 
a perfect fluid solution (at least for the seed metrics considered). Fi- 
nally, starting with Tolman IV, Heint Ha and Durg IV, V solutions, we 
build anisotropic gravastar-like sources with the appropriate boundary 
conditions. 

Keywords: Anisotropic spheres; Tolman solution; Florides solution; Boson 
stars; gravastars. 
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1 Introduction 

In the real world, stars are not made up of perfect fluids. Despite this, it 
seems that perfect fluids are appropriate for describing ordinary astrophys- 
ical objects. However, since of the work of Bower and Liang pQ, anisotropic 
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models have been attracting increasing interest [2] . Firstly, anisotropic fluids 
have been considered in studying objects with hight redshift such as quasars 
[3]. Local anisotropics can arise from modeling very compact objects such as 
neutron stars. Further, anisotropies can lead to important modifications of 
the physical parameters of the fluid [H [5J. In particular, these can affect 
the critical mass and the stability of stars [8] and also the redshift. In 
this context, particularly interesting is the case of boson stars [HUH], where 
naturally anisotropies come in action. These are hypothetical stars, born of 
an idea of Kaup [15], constituting a bosonic particle solution of the Klein- 
Gordon equation coupled with general relativity. Bosons are described with 
a complex scalar field $(r, t),: <3?(r, t) = J(r)e~ lat , with a the frequency and 
(r, t) the usual coordinates in a static spherical space-time. The function 
J(r) is a measure of the anisotropy of the star. These stars always have 
(P± — P r ) < 0, with P± and P r respectively the tangential and the radial 
pressure, and (P± — P r ) vanishes at the centre r = and near the approx- 
imate radius of the star (see and references therein). From the consid- 
erations above, it is evident that there is a need to investigate anisotropies 
in general relativity. In particular, it is useful to build anisotropic metrics 
that differ from a given seed metric by free parameters. In the literature 
[U [UJ HS1 [J7], methods have been developed to build explicit anisotropic 
spheres. These methods require, before solving the Tolman equation [19] . 
an ansatz regarding the expression of the anisotropy (P± — P r ). Generally, 
only solutions with constant density E or with E tx are pratically avail- 
able [U CD] . In other cases [TBI Ej , anisotropy factor is calculated by means 
of an integration involving the chosen seed perfect fluid configuration. Also 
in this way, only particular configurations are available. Recently |18j . all 
the anisotropic static spherical solutions have been obtained via quadratures 
by specifying two generating functions . 

Motivated by these facts, we derive a simple procedure for obtaining work- 
able regular anisotropic metrics with the anisotropy factor depending on 
arbitrary regular functions satisfying some boundary conditions. We point 
out that this paper is not principally devoted to the generating method, but 
rather to the physical consequences of a departure from a physically sensible 
configuration. In particular, we are interested in possible phase transitions 
of the matter caused by a "small" departure from a given starting space- 
time. We show as a "dynamic" describing this deparure can be introduced. 
We present our procedure in detail with two seed space-times: the Florides 
|20j class of solutions with vanishing radial pressure and a particular Tolman 
[2T] perfect fluid solution. Further, we show that the technique can also be 
applied, at least, to other perfect fluid seed metrics than the Tolman IV 
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one, i.e. Heint Ha, Durg IV and Durg V. In section 2 we describe our simple 
procedure together with the regularity conditions on the metric coefficients. 
In section 3 the technique is applied to the Florides solutions and in section 
4 to the Tolman IV solution and to the Heint Ha metric. In section 5, start- 
ing with the Tolman IV and Heint Ha metrics, we build a simple gravastar 
model. Gravastar models can also be built (at least) starting with Durg IV, 
V spacetimes. Section 6 presents some final remarks and conclusions. 



2 General Facts 

Our starting point is the line element appropriate for static spheres: 

ds 2 = A(r)dr 2 + r 2 [d0 2 + sin 2 6>c^ 2 ] - B(r)dt 2 . (1) 

The most general expression for the energy-momentum tensor T ab is: 

T ab = E{r)V a V b + P ± (r)[W a W b + L a L b ] + P r (r)S a S b , (2) 

where E(r) is the energy-density, P±(r) the tangential pressure and P r (r) 
the radial pressure, with (Einstein's equations are G ab = —T ab ) 



V a 



W a = [0,r,0,0}, 
L = [0,0, r sin 0,0], 



(3) 



S a 



y/A{r), 0,0,0 



In terms of A(r),B(r), the Einstein's equations reduce to: 
rA r + A 2 - A 



E = 



r 2 A 2 



1 



r 2 A 



4rA 2 B 2 



, B r 

1 - A + r^- 
B 

X(r), 



X(r) = 2B, r AB - 2A, r B 2 + 2rB^ r AB - rA^ r B^ r B - rB 2 r A, 

1 



P ± - P r = S(r 



Mr), 



4r 2 A 2 B 2 

r(r) = 2rA, r B 2 + 2rB, r AB - 2r 2 B^ r AB + r 2 B^A^ r B - AA 2 B 2 + 
AAB 2 + r 2 AB 2 r , 



(4) 

(5) 
(6) 

(7) 
(8) 
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where subindices with comma denote partial derivative. 

Now, suppose that A(r), B{r) depend smoothly on certain parameters Oj, b{ 

respectively in such a way that, setting a,{ = 6j = 0, we have 

A(r)=A (r) , B(r) = B (r), (9) 

A , B being given seed solutions. 

Under the assumption that the seed solution is regular, satisfies all energy 
conditions and is smoothly matched to the vacuum Schwarzschild solution, 
we show that the solution A(r),B(r) also is, provided that aj,6j is chosen 
appropriately and that, by inspection of (|5])-([8|), the following conditions are 
fulfilled: 

A(r) >0Vr e [0,R] , A(R) = A (R), (10) 
A(0) = 1 + ar n + o(l) , a e R , n > 2, 
B{r) > Vr G [0,R] , B(R) = B (R) , B, r (R) = B , r (R), 
B r (0) = (3r n + o(l) , p e R , n > 1, 

i? being the radius of the star. Further, if the seed metric is a perfect fluid 
solution or anisotropic matter with P ± (R) = 0, then the generated solution 
will have S(R) = 0, provided that 

A tT (R) = A 0}r (R) , B„{R)=B^ r {R). (11) 

Finally, all energy conditions and physical reasonability for usual ordinary 
matter require that 

P r > 0, P± > 0, E > , E - P± > , E - P r > Vr € [0,R]. (12) 

For usual matter, all of the conditions (I10j) - (|12j ) must follow. Conversely, we 
can build solutions with non-usual matter, by dropping some of the condi- 
tions (|12p . In the next two sections, we show that the procedure described 
above can effectively generate physically reasonable anisotropic fluids. 
As seed metrics we take the Florides solution without radial pressure and 
with constant energy-density and the perfect fluid solutions called Tolman 
IV and Heint Ha. 

3 Example 1: Florides seed solution 

As a first example, we take as seed metric the Florides |20| class of solutions 
representing anisotropic static spheres with vanishing radial pressure. The 
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most interesting feature of these solutions is that the tangential pressure 
is sufficient to prevent the gravitational collapse of a dustlike distribution 
of matter. Further, these solutions describe the interior field of Einstein's 
clusters where the tangential pressure is generated by a gas of noncollisional 
particles moving on spheres. From Q we see that these solutions have 

A (r) = l_ = l + r ^, (13) 

uw , 2m (r) Bn 



m (r) 



i r 

- / E (u)u 2 du , m {R) = M, 
2J 



where mo(r) is the gravitational mass inside a radius r. As a title of example, 
we consider the solution in [20] with constant energy-density given by 



» / \ 3 „ 6M . , E 2 r 2 

Mr) = ^~ ^ , E=— , P ±(r)- 



3-r 2 £ ' i?3 ' u ^ v ; 4[3-r 2 £ ] 



3 -1 
?2 77i\ 2 / ^2l^\ 2 



P 0r = , B (r)=(l-^V(l-^l . ill! 



We take 

A = A , B = B e H ^ , E = E . (15) 

The choice for B{r) in (|15p is not restrictive and pemits simpler computa- 
tions. With (|15p . for P r we read 

rvl 3r v ; v 7 

The expression (|16|) implies that, in order to obtain physically reasonable 
solutions, the function H(r), adjoining (fl~0|) - (fT2"|) . must be monotonically 
increasing. 

As an example, we choose 

jp-A > a>0 - 

With (HZD we get 

P r = ^(3-Er 2 )(r 2 -i? 2 ) 2 . (18) 

In this way, P r is regular, positive, vanishing at the boundary r = R and 
monotonically decreasing. Therefore, energy conditions impose E — P r > 0, 
i.e. 

Ml 

»<- = -. (19) 
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Regarding P± , we have 



P± = Po± + ^ (R 2 - r 2 ) (R 2 - 3r 2 ) (R 3 - 2r 2 M) + 



2J2 



9a r 



R 2 - r 2 ) 4 (R 3 - 2r 2 M) . (20) 



The first condition to be imposed for (|20p is that P± > 0. Since the third 
term on the right hand side of (120j) is positive inside the star, the positivity 
of P± can be obtained by requiring that the minimum of the two first terms 
in the right hand side of ([20]) be positive. With a simple algebra, we get 



< -243(3/c 3 + 8k 2 - a)k 6 _ * 

a ~ (6k 3 - 8k 2 + a) 2 (3k 3 + 2k 2 - a)(3k 3 - Wk 2 - a) ~ ° ' 
a = k 2 V$k 2 -24ifc + 28. (21) 



The inequality (|2ip contains (|19p . It is interesting to note that expression 
(|2ip is monotonically decreasing in k for fc > 4, with a — > for — > oo. 
This indicates that anisotropy corrections to the seed metric considered are 
relevant for compact objects. 

Regarding the behaviour of (f20j) . for a < |q, P± is monotonically increasing 
(remember that Pq± is monotonically increasing). Further, for a > j^, Pj_ is 
monotonically decreasing up to a certain value r* where there is an absolute 
minimum and is monotonically increasing for r > r*. 

Thus, for a ^ jq, the term i?(r) drastically changes the behaviour of P±. 
The condition 12 — P± > is satisfied in the range (|2ip . For the anisotropic 
factor S(r) (see (jSj)) we get 

^(^-^'(R^-lrHl). (22) 

For a < jq, expression (|22|) is positive and monotonically increasing inside 
the star. Conversely, when a > j^, we have always 5(0) = but S(r) is 
now negative up to a certain radius r' where S(r') = 0. For r G (r', R], S(r) 
becomes positive and monotonically increasing. Therefore, the physics of 
the system changes for a sufficiently large. Note that with a negative value 
for S(r) we have a "soft"' core, while for S(r) positive we have a "hard" core 
for the star. Further, a negative value for S(r), vanishing at the centre r = 0, 
is typical of boson stars (see [H]). For such stars S(r) = —J 2 r /A(r), where 
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$(r, t) = J(r)e~ tat is a complex scalar field satisfying the Klein-Gordon 
equation and a is the frequency. Obviously, since Klein-Gordon equation 
must be satisfied for boson stars, we do not claim that the solution found is 
an exact solution for such stars, but only that solutions with a behavior for 
S(r) related to that of boson stars can be obtained in our context. Physically, 
the parameter a induces a phase transition for the matter inside the star. 
In fact, when o > jq, there is a region filled with usual matter, i.e. for 
r G (r',R\, but the region r G [0, r'] is filled with unusual matter with 
S(r) < 0. Remember that S(r)/r is a measure of the force exerced by the 
pressure. In the next section we show that a "dynamic" can be introduced 
by means of a time dependence of the parameter a. For the possibility of 
existence of mixed boson-fermion stars, see |22| and references therein. 
Althought we have analyzed only a particular class of generating functions 
H(r), the reasonings above hold for more general choices, provided that 
H[r) isw monotonically increasing, regular and satisfies all the appropriate 
conditions. For example, we can take H(r) = a(r 2 /R 2 - lf + , with a > 
and n a positive integer. The physical features remain the ones depitced 
above, with an expression for a* similar to (|2ip with k = wf^n ■ 



4 Example 2: Tolman seed solution 

As a further example, we take the well know Tolman IV perfect fluid solution 
derived in |21|: 



R 3 (R 3 -3R 2 M + 2r 2 M) 

o(r) ~ (R 3 - 3R 2 M + r 2 M)(i? 3 - r 2 M) ' ( ' 

3M Mr 2 
B (r) = 1- — + 



R R 3 ' 

_ 3M(2R 6 - 9R 5 M - 7R 2 M 2 r 2 + 2r 4 M 2 + 3R 3 r 2 M + 9R 4 M 2 ) 
° ~ R 3 {R 3 -3R 2 M + 2r 2 M) 2 

3M 2 {R 2 -r 2 ) 



R 3 (R 3 - 3R 2 M + 2r 2 M)' 



Regularity and energy conditions follow for (|23p . provided that R > 3M. 
Note that, since -Eq is monotonically decreasing, we can express r 2 in terms 
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of Eq. After some little algebra, we obtain 

r 2 = F(E ), (24) 
12M 2 R 5 E + 9R 3 M 2 - 4E R 6 M - 21R 2 M 3 + R 2 f3 
( o) ~ 2(4E R 3 M 2 - 6M 3 ) ' 

= V3M 3 (i? - M)(69M 2 - A8ME R 3 - 21RM + 16E R 4 ). 

Putting F(Eq) in the expression for Pq, we get Po = Po(Eo)- 

First of all, we examine the solution A = Aq, B = Boe H ( r \ It is interesting 

to note that the related expression for P r reads 

Pr = Po + ^- (25) 

Therefore, starting with a know perfect fluid metric and by setting H )T = 
— ArPo, we obtain solutions with vanishing radial pressure. 
As an initial imput, we again take H = a(r 2 /R 2 - 1) , a > 0. With this 
choice, we get 

Pr = Po(r) + -^(R 2 - r 2 ) 2 (R 3 - r 2 M) (R 3 - 3R 2 M + r 2 M) , 
Q = R 3 - 3R 2 M + 2r 2 M. (26) 

For R > 3M and Va > 0, P r is monotonically decreasing. Thus E — P r > 
for 

(2RM + 13^ - 11M 2 ) , 
~ 2{R 2 + 3M 2 - ARM) ' 1 1 

For Pj_ we obtain 

_ 3 7 (r)q (it! 2 - r 2 ) 
P± - Po + _ R9 Q 2 (28) 



i? 2 - r 2 ) 4 (R 3 - r 2 M) (Mr 2 - 3R 2 M + R 3 



9r 2 a 2 
R^Q [ 

7(r) = A0R 8 r 2 M - 12MR 10 + 18r 8 M 3 - 6R 9 r 2 + 93R A r A m 3 + 
18M 2 R 9 + 6R 3 M 2 r 6 + 13i? 5 r 4 M 2 - 76R 2 M 3 r 6 - 57R 7 M 2 r 2 + 
2R U - 12R 6 Mr A - 27M 3 R 6 r 2 . 



By inspection of (j28l) . we see that 7 is monotonically decreasing and gets 
its mininum (negative) value at r = R. Further, the third term on the right 
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hand side of (|28|) is always positive inside the star. As a result, the values 
of a for which Pi > are 



M 2 R e (R s -3R 2 M) _ (fe-3) _ 

°- — — ~ A { k-m-if- a - 



The condition (|29[) contains the condition ()27|) . The expression (|29p is a 
monotonically decreasing function of k = R/M for k > ■ As we shall 

see later, stability requires k > 4.661. Consequently, according to [2J, the 
role of anisotropies becomes more relevant for stable compact objects. 
Further, E — P± > and P± is monotonically decreasing in the range (|29p . 
For S(r) we have: 

P ± - P r = S(r) = ^r 2 A (R 2 - r 2 ) + (30) 

^^{R 2 - r 2 ) 4 {R 3 - 3R 2 M + Mr 2 ) (R 3 - Mr 2 ) 

A = -4R 9 + 4r 4 M 2 i? 3 + 57r 2 M 3 R 4 + 14M 3 r 6 - 8r 2 MR 6 + 
24Mi? 8 - 54r 4 M 3 R 2 + 9r 2 M 2 R 5 - 33M 2 R 7 - 9M 3 R 6 . 



By a Taylor expansion of (|30p near r = we get 

_ 3ar 2 [3 J Ra(fi-3M)-4P 2 + 3M 2 ] , 
i(r) " &(R-3M) + ( '■ [ ' 

From ([BT|) we see that S(r) > near r = for a > ^^zfy- Thus, thanks 
to ||22|l, we have S(r) < 0,Vr G [0,fl] with 5(0) = S"(i?) = 0. Since 5(r) 
is negative, the star has a "soft" core. The behaviour above for S(r) is 
typical of the (hypothetical!) boson stars [9l [14]. To best represent such 
objects we could require that (remember that a boson star has no well- 
defined radius [11]) S(r) vanishes exponentially outside the star reading a 
maximum negative value near r = R. This can be obtained, for example, 
by dropping condition (fTT|) and by posing 

H(r) = -a(r 2 /R 2 - 1)V^ , n > 2 , a > 0. (32) 

However, the behaviour of P r ,P± is very similar to the one depicted above. 
Concerning the equation of state of our generated solution (j26"l) - ([28j) . thanks 
to ([Ml), we have P r = P r {E),P ± = P±(E), with P r>E £ [0, 1),P±, E e [0, 1), 
where 

~6M(R-2M) 3M{2R-3M)~ 



E £ 



R 3 (R-M) ' R 3 (R-3M) 



(33) 
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since E is monotonically decreasing. 

For example, for k = 10, a = a* , (Pr,E) max — 0.33. Therefore, the sound 
velocity of the fluid remains always less than 1. For the seed solution, 
(Po,E) max — 0.24. However, note that for Tolman IV solution, Pq^e is a 
monotonically decreasing function of E. Conversely, with respect to the 
generated solution, P T) e is a monotonically decreasing function of E for (ap- 
proximately) a G [0, ^j-), while is a monotonically increasing function for 
a G [%-,a*]. Note that, since E is a monotonically decreasing function of r, 
we expect, on general grounds, that the sound velocity decreases for increas- 
ing r. This desiderable behaviour is obtained in our solution in the regime 
where anisotropies are most relevant. Therefore, for a G [^-,a*] our gener- 
ated solutions satisfy all the fundamental features of physicity requested in 
the literature. 
For the redshift we have 

_H 

e 2 

1. (34) 



With the starting choice for H and for small a, we read 

, = ,o + ^=(l-^) 3 + o(a). (35) 
The maximum redshift for our solution is attained at the centre: 



1 + - 

2 -l + o(o). (36) 

1 _ 3A£ 
1 R 



If we impose the maximum redshift possible for stable boson stars [14] . i.e. 
Zmax — 0.687, for the critical value k cr u we get 

3 

k > k cr n , k cr n — - , . o • (37) 

1 -0.351(1 + a/2) 2 

From the reasonings belove equation (I29p . the maximum value of a* is at- 
tained at k = ■ By putting this value in (|37p , we obtain the estimation 
k cr n — 4.661, and thus k > k C rit f° r stability. Note that the redshift is bigger 
than that of the seed Tolman IV solution. 

By taking B(r) — ► B(r)T(t), the field equations do not involve derivatives 
with respect to the time coordinate. As a result, if we take a — > a(t), all the 
expressions given above remain unchanged. Consequently, we can add "dy- 
namic" to our generated solutions by choosing a generic time dependence 
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for a, provided that a(t) < a*. For example, we can take 

a*t 2 
o(t) = — 



T 



1 



(38) 



r being a constant. Physically, the term a(t) could be interpreted as the 
"ignition" of some scalar field, represented by the function H(r,t), at t = 0. 
This scalar field can generate a phase transition ending with the transforma- 
tion of the starting "Tolman" star in a boson- like star. The considerations 
above suggest the possibility that boson stars can form by a suitable "per- 
turbation" of a perfect fluid (at least for the seed metric considered). For a 
review for this problem see |24} [25] . Finally, also for the solution depicted 
in section 3, we can introduce a time dependence for a given by f)38[) with 
a* given by (f2~Tl) . 

We studied the physical properties of the generated solution in terms of 
the generating functions H{r) satisfying conditions (I10j) -(|12p. The results 
obtained cannot be summarized as a theorem, but, almost, only necessary 
conditions can be done. For example, from a lot of functions H(r) consid- 
ered it seems that, in order to have S(r) < 0, H(r) must be monotonically 
increasing. Conversely, in order to have 5(r) > 0, H{r) has to be chosen 
monotonically decreasing. 

For example, for H(r) = a(r 2 /R 2 - lf n+3 ,H(r) = -a{r 2 /R 2 - l) 2n+4 ,a > 
0, where n is a positive integer, we have 5(r) < inside the star, while the 
same expression with a negative leads to 5(r) > 0. But monotony is a neces- 
sary but not sufficient condition. For example, for H(r) = a(r 6 /i? 6 — 1) , 
5(r) changes sign inside the star. Generally, if H(r) is not monotonic inside 
the star, 5(r) changes sign. For example, if H = ar 2 /R 2 (r 2 /R 2 — 1) , a > 
0, then S(r) > up to a certain r = r* ("hard" core), while S(r) < 
("soft" core) for r G (r*,R] and 5(0) = S(r*) = S(R) = Conversely, if 
H = —ar 2 /R 2 {r 2 /R 2 - l) 2n+3 ,a > 0, then S(r) < up to a certain r = r*, 
while S(r) > for r G (r*, R] and 5(0) = S(r*) = S(R) = 0. 
We now study the generated solution with B = Bo, A = Aqc g (without loss 
of generality). For example, we could take G = br 2 /R 2 {r 2 /R? — l) 4 ,fr > 0. 
Also in this case, we can find a limiting value for 6* such that for b < b* 
the generated solutions satisfy all the conditions (fTU|) - (fT2j) with E(r),P r 
monotonically decreasing. Further, for b < \, P± is also monotonically 
decreasing. For example, if k = 10, we found b* ~ Further, for all the 
values allowed for b, S(r) > and 5(0) = S(R) = 0. If we put, in the exam- 
ple above, b — > — b with b > 0, then we read 5(r) < and 5(0) = S(R) = 
with the same range of variability for b. 
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Furthermore, we have studied the generated solutions by considering a lot 
of G(r) satisfying conditions (fT0l) - (fT2|) , Also in this case, necessary and suf- 
ficient conditions cannot be imposed. However, it seems that a necessary 
condition to have S(r) < is that G(r) < 0, while a necessary condition to 
have S{r) > is that G(r) > 0. 

To conclude our study, we consider A = Aoe G , B = B$e H . 
The results obtained are: 

• If H(r) is such that S(r) < for G(r) = and G{r) is such that S(r) > 
for H{r) = 0, then for a ~ a* and b < b* /n for a certain n > 1 de- 
pending on the explicit expression chosen for (H,G), H{r) dominates and 
S{r) < Vr G [0,12] and S{R) = S(0) = 0. Conversely, if a < a*/n and 
b ~ 6*, then G(r) dominates and therefore S(r) > Vr G [0, J?]. 

• If i?(r) is such that S(r) > for G(r) = and G(r) is such that 
S(r) < for i?(r) = 0, then for a ~ a*, b < h*/n, H(r) dominates and 
S(r) > Vr G [0, R]. Conversely, if a < a*/n and 6 ~ 6*, then G(r) domi- 
nates and S(r) < Vr <G [0, 12]. 

• If H(r) is such that 5(r) > for G(r) = and G{r) is such that S{r) > 
for H(r) = 0, then 5(r) > Vr G [0, 12]. 

• If H(r) is such that S(r) < for G(r) = and G(r) is such that 5(r) < 
for H{r) = 0, then S(r) < Vr G [0, 12]. 

• In all the other cases, S(r) changes sign inside the star. 

Obviously, the values of a*,b* depend on the chosen generating functions 
(H, G) , but generally we get expressions similar to (|29p . 
From the discussion above, it follows that, in particular, configurations with 
S(r) < can be generated also by varying the energy density of the seed 
perfect fluid solution. 

The results obtained with Tolman IV, in relation to the possibility of build- 
ing interior solutions with S < 0, do not represent a "lucky" case. The same 
"phenomenon" can be easily obtained with other perfect fluid seed metrics. 
As a further example, we can take the seed Heint Ha metric ( see |23|). The 
metric functions are (12 > 3M): 



B = K(l + ar 




3ar 2 ( 1 + C(l + 4ar 2 ) 2 



) 



M 



K 



2 I 1 + ar 2 

(312 - 7M) 3 



a = 



R 2 (3R-7M) 



2712(12 - 2M) 2 ' 




(39) 
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With the metric ([39]), we have S(r) < Vr G [0, i?] with i?(r) given by (HU) 
(with a > 0) and 5(0) = 5(i?) = 0, privided that: 



24k 2 - 65k + x /3g=H(9ife - 24) 

• < ' (40) 



By using the less restrictive condition (|29p . it is follows that, for a > 0, in 
adjoint, our generated solutions satisfy all energy conditions and Vr G [0, i?] 
are regular and E(r), P r (r), P± (r) are positive and monotonically decreasing. 
Thus the seed metric Heint Ha can also be used to build models recalling 
boson stars. Further, the procedure above for building boson-like stars can 
be successfully applied to Durg IV and Durg V seed metrics with H (r) given, 
for example, by (fT7|) . (see[23j). 

5 Example three: Building gravastars solutions 

Gravastars (see \26\ I27j) represent an alternative to black holes. These are 
described in terms of an interior with the equation of state E = —P (de 
Sitter) matched with a thin shell with E = P. In [28] a more generalized 
class of gravastars is discussed for which no thin shell is needed. The matter 
is composed by an anisotropic fluid with a continuous expression for P r with 
the following boundary conditions (see |28|): 

1) Vr G [0,R], E>0 

2) E(0) = -P r (0) = -P±(0) (only at the centre), P r , r (0) = 0. 

3) For r* ~ 2M, there exists a positive maximum for P r . For r G [r* ,R), P r 
is monotonically decreasing and the stars is filled with "ordinary" matter. 
For r G [0,r*), P r is monotonically increasing. Further, S > 0. 

4) There exist exactly two points at which P r = 0. 

All these conditions can be easily fulfilled. As an example, we consider again 
the Tolman IV solution together with 

(2 \ ^ 
l- r v A , a>0, 6G [0,R). (41) 

In this way, the solution for r G [0, b) is composed of anisotropic matter, 
while for r G [b, R] the matter content is that given by Tolman IV solution 
(|23p and the two metrics are smoothly joined at r = b. Condition 1) is 
obviously satisfied. Conditions 2) follow by setting: 

_ Mb\R-M) 

~ R\R - 3M) " 1 ' 
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Regarding conditions 3), we must to chose b ~ 2M. It follows that S > 
and 5(0) = S(b) = 0. Finally, the behaviour for P r is the one required by the 
conditions 3), 4). For example, by setting b = 2M, R = 10M, we found r* = 
1.9986M < b. The graph of P± is similar to the one for P r and for P r and P± 
follow expressions similar to those given by ([25]) . ([25]) . Obviously, at r = R 
our solution is smoothly matched to the vacuum Schwarzschild solution. 
With the simple procedure described in this section, specific workable models 
for gravastars can be built simply by setting an appropriate perfect fluid 
solution in the region r G [6, R] . 

We can use other seed metrics than Tolman IV. As a further example, we 
can use Heint Ha metric given by (139p . For this seed metric, it is easy to 
see that all the conditions 1 — 4 can be fulfilled with ([3T]) . provided that, 
instead of (1321) . we take: 



Mb 2 



3R + 



|^( 3j R-8M) 



2R 3 (3R-7M) 



(43) 



with b ~ 2M. We do not enter in the details, but only mention the fact that 
the procedure above to build gravastar-like models can also be obtained, for 
example, with the Durg IV and Durg V seed metrics, with respectively 



4M6 2 



21i? 3 (4i?-9M) 2 

5Mb 2 
336i? 3 (5i?- 11M) ; 



4R-9M 

5{R-M) 
5R-11M 



117 MR 



(44) 



+ 29161M 2 i? - 265WMR 2 + 6025i? 3 



224M 2 - YJ2MR + 32i? 2 , 



s = 2375iT + 52425M 2 i? - 44800M 3 - 19800Mi? 2 



With the relations (|44p . also for Durg IV, V all the conditions 1 — 4 can be 
sasfied together with b ~ 2M. 



6 Conclusions 

Motivated by the increasing interest in anisotropic fluids, we presented a 
simple procedure for building physically reasonable spacetimes with local 
anisotropies starting from a known regular physically sensible seed metric. 
The principal aim of this paper has been the construction of a simple algo- 
rithm for exploring the physical consequence derived by a departure from 
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a given configuration by means of given parameters. We studied in detail 
two physically interesting cases by taking as seed spacetimes the Florides 
solution without radial pressure and with constant energy-density and the 
Tolman IV perfect fluid solution. In particular, we showed that spacetimes 
with a behaviour of the anisotropy factor S(r) closely related to that of bo- 
son stars can be generated starting from a suitable perfect fluid solution (at 
least with the seed solutions considered!). A time dependence can be intro- 
duced in the B(r) metric coefficient simply by taking a — ► a(t). The same 
is not possible for A{r) because, in this case, G r t / and a heat flow term 
appears. Physically, the term a(t) could be interpreted as the "ignition" at 
t = of some scalar field represented by the term H{r) in the metric compo- 
nent gu- Further, starting with an appropriate seed perfect fluid metric, also 
a gravastar-like model can be built with the desired boundary conditions. 
In particular, we have obtained gravastar-like models and boson-like stars 
by using the following seed metrcics: Tolman IV, Heint Ha, Durg IV, Durg 
V. 

Finally, we point out that the technique used in this paper can also be 
used by taking a generic starting metric representing a regular physically 
reasonable perfect fluid with an oculate choice for the functions H(r), G(r). 
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